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1 Introduction 



The complicated vacuum structure of quantum field theories involving nonabelian 
gauge symmetries has been a subject of considerable interest since its consequences 
have been discussed for the first time |[ . In such theories, besides the usual quan- 
tum states built perturbatively on the topologically distinct vacua (characterized by 
a topological quantity, the winding number), there are nonperturbative phenomena 
relating these separate perturbative sectors In particular, there are transitions 
between the states built on neighbouring vacua which are dominated by different 
physical processes, depending on the temperature of the system. These processes 
can be described semiclassically by pseudoparticles, i.e. classical configurations with 
particle-like properties. The most important example for these effects is the elec- 
troweak model in which winding number transitions are related to baryon/lepton 
number nonconservation Most calculations, however, are performed in sim- 

pler, more manageable "toy models" . 

This article presents the evaluation of such a transition rate in a special tempera- 
ture region for a particular but typical model which has the advantage of permitting 
explicit calculation in a nontrivial case. In this introduction, we therefore first present 
the general physical background, i.e. the physics of transitions between topologically 
distinct sectors of a quantum field theory in different temperature regions, and then 
introduce the particular model which is motivated by its special temperature be- 
haviour, but also by some special technical properties. 

1.1 Transition rates due to tunneling or thermal fluctuations 

Depending on the temperature of the system, transition processes between neigh- 
bouring perturbative sectors of a quantum field theory are dominated by different 
physical processes: 

At zero temperature T = 0, the relevant process is tunneling through the energy 
barrier separating the neighbouring vacua. In a semiclassical approximation, this 
tunneling can be described by vacuum instantons 0, 0. Vacuum instantons are 
pseudoparticles which are stable classical solutions of the Euclidean Euler-Lagrange 
equations of the theory with vanishing Euclidean energy as boundary condition, 
^inst = 0. A semiclassical expansion of the Euclidean path integral around the 
instanton yields the rate of tunneling through the barrier separating the topologically 
distinct vacua between which the instanton interpolates as a funtion of imaginary 
time, used here to describe tunneling. The tree approximation of the instanton 
tunneling transition rate is thus exp{— 2150} where 5*0 is the Euclidean action of the 
vacuum instanton ||^. 

At finite temperatures T > 0, there are thermally excited states in the neigh- 
bouring sectors, and starting from some lower temperature limit, transitions due 
to thermal fluctuations over the separating barrier become relevant. These transi- 
tions are described semiclassically by pseudoparticles called sphalerons. These are 
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classical solutions of the static Euler-Lagrange equations with finite Euclidean en- 
ergy £sph which are unstable, i.e. the static Gaussian fiuctuation operator around 



these solutions has one negative eigenvalue U- |T0|- Sphalerons can be visualized 
as "sitting on top" of the energy barrier separating neighbouring vacua. The corre- 
sponding semiclassical transition rate is suppressed by the classical Boltzmann factor, 
exp{— £^sp/i/?t} which one obtains as the tree approximation of the partition function 
path integral which again involves the use of complex time, now in the convenient for- 



malism of quantum field theory at finite temperature [11, f3l. Comparing the leading 



exponential factors of the semiclassical expansion around instantons and sphalerons, 
respectively, one finds that thermal fiuctuations dominate vacuum tunneling in the 
transition rate for temperatures T > Tq =: y^^^- 

A third transition process which has been discussed intensively in recent literature 



14, 15 1 is tunneling from (thermally) excited states, described by a Boltzmann average 



over nonvacuum instantons which are classical solutions of the Euclidean Euler- 
Lagrange equations with Euclidean energy Einst > as boundary condition |16|, p!7| . 
These solutions are periodic in the imaginary time coordinate with energy dependent 
period f = f{Sinst) (which is why they are usually called "periodic instantons"). In 
order to apply the Boltzmann average over the periodic instantons, one has to satisfy 
the saddle point condition f{Sinst) = Pt, so that only periodic instantons with period 
equal to the inverse temperature dominate the nonvacuum tunneling transition rate. 
The Euclidean action of the periodic instanton can thus be written SinstiPr), and 
the usual relation Emst = '^'"^^'^^^ holds. The two different uses of complex time (in 
tunneling and in thermodynamics) thus merge together in the case of the periodic 
instanton. 

For Sinst £sph the imaginary time dependence vanishes and the periodic instan- 
ton reduces to the sphaleron (being time periodic with any period, i.e. static) with 
Euclidean action Sgphif^r) = ^sphPr- The second limit Sinst which should reduce 
the periodic instanton to the vacuum one shows that there are two different types of 
periodic instantons |T^: Those with f (Sinst) oo where one has localised vac- 

uum instantons with infinite period and those with f (Sinst) where localised 

vacuum instantons do not exist. Nonetheless the well-known method of constrained 
instantons [|r^, [T^ (which may be interpreted as the zero period limit of this second 
type of periodic instantons) allows vacuum tunneling even in this case. In either case 

the dilute gas approximation works and one has Sinstd^r) 250. 

It is now an interesting and important question to ask how the temperature do- 
mains in which the different physical processes (vacuum tunneling, thermal tunnel- 
ing, thermal fluctuations) dominate the transition rate, can be continued into each 
other. The crucial point here is the crossover from predominantly thermal fluctua- 
tions to predominantly (vacuum or thermally excited) tunneling. Near the sphaleron 
energy, Ssph — Sinst ^ Sinst, the periodic instanton can be written as a sum of the 
sphaleron and oscillations in its eigenmode with negative eigenvalue uj^. In the 
limit Sinst — > Ssph, the period of oscillation approaches fsph '■= zf' ~- kBT(-) ' 
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Figure 1: Type I theory (the solution that dominates the transition rate is indicated 
by the thick line) 



periodic instanton thereby merging into the sphaleron which has Euclidean action 
Ssphitsph) = slpl at this lower temperature bound. Depending on the model, the 
period f{£inst) for Einst < £sph may be smaller or larger than fgph- 

This fact, together with the two types of low energy behaviour of the periodic 
instantons, results in two different simple types of theories with continuous crossover 
between the temperature domains with different dominant processes |jl5|. In type 

I theories, we have localised instantons, f{£inst) ^"^^^^ oo, the period f is always 
greater than fspt, thus = %f < 0, and 250 > Si~l Then Tq < T^-\ 
thermal fluctuations dominate for T > T'^~^ and nonvacuum tunneling from thermally 
excited states dominates for T < T^~\ and is suppressed by exp{—Sinst{PT)}- The 
crossover between the two regions is smooth. This situation is shown in Fig. |l|. Type 

II theories have nonlocalised vacuum instantons, f{Sinst) 0, the period f is 
always smaller than fsph, thus = > 0, and 25*0 < S'^. This yields 

< Thermal fluctuations dominate for T > Tq, vacuum tunneling for T < Tq 
and the crossover is sharp as shown in Fig. 0. 

Well known examples are the 1 + 1 dimensional Higgs model |jl6|, |19| for a type I 
theory and the reduced nonlinear 0{3)-a model pO| , |14| for a type II theory. Both are 
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Figure 2: Type II theory (the solution that dominates the transition rate is indicated 
by the thick line) 



prominent toy models for the winding number transition in the electroweak model. 
But the electroweak theory seems to be neither of type I nor of type II: It has 
no localised instantons which disagrees with type I properties, and one estimates 
25*0 > SsphTsph for the electroweak constrained instantons, so that the theory is not 
of type II. 

To understand this special property of the electroweak model which is not shared 
by the usual toy models, it is interesting to study simple low dimensional models 
which also can not be classified under the two types above, although such models 
may not be good toy models for the electroweak theory from other points of view. 
One of these models of recent interest is the reduced nonlinear 0{3)-a model in 



1 + 1 dimensions with an added Skyrme term PT|. 

1.2 The reduced 0(3)— cr model with Skyrme term and its 
motivation 

The model we deal with here is defined by the Lagrangian density 

o Z 
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(2) 



in Minkowski spacetime (x , x ) = {t,x) with fields = {(j)""), a = 1,2,3 constrained 
by (p"-(j)°' = 1 (we do not add a Lagrange multipher since we will always work with 
parametrisations respecting the constraint). This Lagrangian differs from that of the 
usual nonlinear a model, C = ^{d^(p"'d'^(p"') ||2^ (comparing the kinetic terms of 
the two Lagrangians, the two couplings p- = Ai have to be identified) with 0"0" = 1 
by two additional terms. The explicit 0(3) — ^ 0{2) symmetry breaking potential 
(which we choose to be V{<p^) = 1 — 0^), being necessary to produce a sphaleron in 
the theory (there is no mechanism of spontaneous symmetry breaking for a models), 
and the Skyrme term with coupling i/"^. Setting z/^ = reduces this model to the 
well-known reduced nonlinear 0{3)-a model which provides good model predictions 
for the electroweak theory pOl- The latter is the original Skyrme-less model which 



has no localised instantons and is a type II theory 



Adding a Skyrme term changes this last aspect completely, 
localised vacuum instanton ifinst in the Skyrmed theory 
numerically ||2T|. Its action is |jl5 



First, there is a 
which can be verified 



S'o = 47rAi + C(AiZ/^). 



(3) 



In this model the sphaleron solution Lp^ph is not affected by the Skyrme term 
because it solves the static Euler-Lagrange equations and the Skyrme term is purely 
dynamical. But the fluctuation spectrum around the sphaleron is modifled and there- 
fore also the one-loop quantum correction in the semiclassical transition rate. More- 
over, the temperature range of validity of this calculation changes. In particular, the 
temperatures To and T^~^ will depend on u. 

It is the purpose of this investigation to calculate the thermal transition rate 
in a semiclassical one-loop expansion around the sphaleron in this Skyrmed theory. 
This involves the evaluation of fluctuation determinants, and one main motivation 
of this paper (besides knowing that the Skyrmed reduced nonlinear 0(3)-cr model is 
a theory neither of type I nor of type II) is that these evaluations can be performed 
explicitly, and one can flnally see clearly the role of the Skyrme term in the thermal 
transition rate. 

Moreover, the technical aspects of the calculational difficulties which the Skyrme 
term gives rise to are presented in detail. These techniques should be of interest 
also in the context of more physically motivated models like the three dimensional 
a model which is motivated by QCD. In the case of the latter, the hope is that the 
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loop contributions could give rise to Skyrme-like terms in the action which would 
stabilize the Skyrmion . 

The layout of the paper is as follows: In Section ^ the notion of a sphaleron is 
reviewed. The NCL technique is discussed in the context of the Skyrmed 0{3)-a 
model, the sphaleron solution is derived and its instability is discussed in terms of 
the static Gaussian fluctuation operators. 

The embedding of the sphaleron of the reduced 0{3)-a model into the full Eu- 
clidean theory is described in Section |^ where the thermodynamical formalism yield- 
ing the transition rate is reviewed, using the path integral version of quantum statis- 
tics [^. What is quantized are the fluctuations parallel and perpendicular to the 
vacuum and the sphaleron, respectively. Since these fluctuations are dynamical quan- 
tities, they may involve the (kinetic) Skyrme term. The effect of the latter on the 
model is to produce an extra measure factor which modifies the scalar product in 
the Hilbert space of the perpendicular fluctuations around the sphaleron. Therefore 
the evaluation of the path integrals is also considered in some detail in this Section, 
taking care of the path integral measures. 

Sections ^ and ^ deal with the fluctuations around the vacuum and the parallel 
fluctuations around the sphaleron which are already well known since there is no 
influence of the Skyrme term on this part of the theory . 



The new effects of the Skyrme measure factor on the perpendicular fluctuations 
around the sphaleron are investigated in Section 

In Section the results of the calculations are summarized to yield the ther- 
mally activated transition rate. The Skyrme-less limit is discussed, as well as the 
temperature range of validity of the calculation. 

The physical imphcations of this analysis are then discussed in the concluding 
Section ||. 

Finally, some mathematical steps and calculational details are given in the Ap- 
pendices. 



2 Saddle— point configurations of the static energy 
functional 

The solutions we consider are static, cj) = (j){x), and result as extrema of the static 
energy functional 

S[$] = J Hstatdx (4) 

where Tigtat is the static part of the Hamiltonian, defined as the 00-component Ti. = 
Too of the energy-momentum tensor T^,^ := 11^0" — g^u^, := 0(0^, or obtained 
alternatively from the well-known formula 

n:=u^-c, n:=^ (5) 
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Although one cannot write down the Hamiltonian in terms of fields and conjugate 
momenta explicitly, it is easy to see that 



H 



Ai 



stat 



(6) 



For the special potential V{(1)^) = 1 — 0^, the vacuum configuration of £[(j)\ is thus 
given by (^0 = (0,0, 1). 

We parametrize the fields so that they satisfy the constraint (jp' = 1. There are 
at least three appropriate parametrizations of the static fields. 



(t){x) 



/ smf^{x) 
\ cosf^{x) 



(j){x) 



/ sin (x) COS g"^ (x) 
sin g^{x) sin g^{x) 
\ cosg^{x) 

^ sin {x) sin h"^ (x) \ 

sin h^{x) cos h^{x){cos h?{x) — 1) 
y sin^ cos /),^(a;) + cos^ j 



(7) 



(9) 



In terms of the new parameter fields /(x), g{x) and /i(z), the static energy functional 
(I) is £[0] = £f[f\ = £g[g] = £h[h] with 



m 

£h[h] 



X,{f'ix),Fifix))f'ix)) + XoV 



cos f^{x) 



dx 



Xi{g{x),G{g{x))g'{x)) + XqV (cosg\x 



dx 



Xi{h\x),H{h{x))h'{x)) + XoV (sin^ h\x) cos h'^{x) 



cos^ h^(x 



with the matrices 
F{f) = 

G{g) = 
Hih) = 





1 



l + (/2)2 

(i+ITW / 
1 \ 
sin^g^ J 

cos^ sin^ /i^ + (1 — cos /i^)^ sin cos sin h? 
sin h} cos sin sin^ h} 



(10) 

(11) 

dx 
(12) 

(13) 
(14) 

(15) 



The /i-parametrization is usually used to find the saddle-point of the static energy 
functional (^ with the special potential V{4>^) = 1 — 0^. Treating h^{x) as a constant 
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parameter h^{x) = t] and writing h'^[x) = f{x), the parametrization yields a 
noncontractible loop (NCL) |26[| 



^ smr]smf{x) 

sinr^cosr^(cos/(x) — 1) 
\ sin^ ?7 cos /(x) + cos^ ?7 



(16) 



in the space of static field configurations on which S is defined. The NCL starts and 
ends in the vacuum. The energy of the configurations along this loop is given by 



Ai sin^ 7] 



-(/'(a:)r + A.^(l-cos/(x)) 



dx 



(17) 



which has a maximum in the r^-direction for = |. Minimizing the resulting energy 
function by variation of f{x), 



Ai 



Sf{y) 



ifix)y + fi\l-cosfix)) 



dx = 



yields the well-known Sine-Gordon equation 

f"{x) - ii^smf{x) = 0, 



/i 



Ao 
Ai' 



(19) 



Assuming as boundary condition that for \x\ oo, 0(x) approaches the vacuum 0q, 
i. e. /(x) — ?^ 0,27r, the solution of (|19D is 

fix) = 4 arctan (^e^'^^^"^'")) . (20) 

From the NCL technique [^], we expect the classical solution which in this 
parametrization is given by = |, f{x), i.e. 0sph{x) = 0i{x) = (sin /(x), 0, cos /(x)), 
to be a saddle-point of the energy functional £[<j)]. This can be verified by consider- 
ing the second variational derivative of S[<p] which should be taken in a convenient 
parametrization of (p to avoid separate treatment of the (f)'^ = 1 constraint. Using 
the parametrizations (^,|^, ^ we thus have to analyse the eigenvalues of the second 
variation of Sf[g\, Sg[f] or Sh[h] taken at the classical solution <^i(x), which means 
/i(x) = fix), f\x) = 0; g\x) = /(x), g\x) = or h\x) = f, h\x) = f{x) 
respectively: 



5r\x)6f^{x')_ 
5g^{x)5gi{x') 
5h'{x)5W{x') 



ifjj^x') dx' = 00 ip^f{x) 



tp^Jx)dx' = u ipli^) 



iljl{x)dx' = u ^l{x) 



(21) 
(22) 
(23) 
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with i,j = 1,2. These eigenvalue equations are equivalent. The /-parametrization 
yields directly decoupled equations in a simple form. Eqs. (|23|) have to be decoupled, 
whereas (p^ ) are decoupled, but can be simplified by a substitution. This is why the 
/-parametrization is preferred in the following, where we will modify our notations 
slightly to f = (/'•"\ Z^"'"'') to denote the two parameter fields. We will see in the 

next section that describes fiuctuations parallel and f^-^^ those perpendicular to 
the sphaleron Lpi{x). In this new notation, each of the eigenvalue equations (PT|-pBD 
may be written as a system of two decoupled static fiuctuation equations, i.e. 

The second equation (^) has one negative eigenvalue (w(^)j = -3^2. ^hus the 
static solution Lpi (x) is always unstable in the space of static field configurations. It 
is important to note that we never left this space. The sphaleron is a static object, 
and one can not decide whether it is Minkowskian or Euclidean. 

The notion of a "static fiuctuation operator" refiects the fact that one can derive 
the stability operator ( [lO|) (which is usually called the "static Gaussian fiuctuation 
operator") by expanding Sf[f] in small "fiuctuations" (which are static in this con- 
text!) around the classical solutions, and similary for the other parametrizations. 
These fiuctuations become dynamical quantities if one embeds the sphaleron into 
the full, time-dependent theory, either in its Minkowskian version (|ip or in the Eu- 
clidean version which one obtains from the former by a Wick rotation in the time 
coordinate, t ^-^ r = —it. As the sphaleron itself does not depend on time even after 
the embedding into the full, e.g. Euclidean theory, one may impose periodic bound- 
ary conditions in the imaginary time coordinate on the fields. Then one can treat 
the result as a field theory at finite temperature which is defined via the imaginary 
time period Pt = ('^^ set fee = 1), and the usual path integral formalism of 
quantum statistics is applicable, including the calculation of thermally activated 
transitions over energy barriers. 



V'^''^^), ^11 = -^ + /^'{^"2sech2(^x)} (24) 
^Y^ij^^\x), = -— + /i' {l - 6sech2 (/ix)} . (25) 



3 Finite Temperature: Sphalerons in Euclidean 
spacetime 



Prom the general theory of thermal transition rates [^, ^ , we know that transitions 
due to fiuctuations at finite temperature may be decribed by the imaginary part of 
the free energy F, 

r = ^-^ImiPrF) (26) 

TT 
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where g is a. damping constant. As usual, F is related to the partition function Z by 
Z = e"^"^^, so that the partition function of the system at temperature T has to be 
computed, using the path integral version of quantum statistics [EH] [O], i.e. 



i=l 



i=l 



1 e 



(27) 



Here we start from the Hamiltonian version of the partition function path integral 
(we emphasize this also by the use of the symbol [d . . .] for the path measures): 
To obtain the Lagrangian version, one has to integrate out the momenta in the 
Euclidean action of the system. 



dr / dx i-iU(f) + n{(j),(j)',Il) 

J-oo \ J 



{2i 



where Ti is the Hamiltonian density (^) of the system. 

If the kinetic term of a field theory Lagrangian is quadratic in the time derivatives, 
e. g. for a field w{t,x) 



C{w,w',w) = -M{wY —U{w,w') 



(29) 



(which means n = Mw), the momentum integration (which in this case is a standard 
Gaussian integral) yields the well-known Feynman formula pSf 



Z = J V{w}e 



Jq^ d,''' J"^ dx Ce(w,w' ,w) 



(30) 



(P{. . .} denoting Lagrangian path measures which depend on the mass M) with the 
Euclidean Lagrange density 



Ce{w, w', ti)) = \m{wY + U{w, w') 



(31) 



where dots now denote derivatives with respect to r. Of course one can also compute 
the Euclidean Lagrangian of our Skyrme model from its Minkowskian version (ffl). 



l\2 



1 

+ 2 



(32) 



but we can not simply use it in formula (^) since the ^-dependence in (|l|) is quadratic 
but involves a field-dependent mass factor, so that the momentum integration has 



to be analysed carefully as it will yield a measure factor |T2|. 

Here we perform a perturbative evaluation of Z around two extrema (there may 
certainly be more contributions): Z ^ Zq + Zi. Zq is evaluated around the vacuum 
00 = (0,0,1) and thus contains the usual perturbative vacuum fluctuations. We 
know that this part has to be real. Zi is evaluated around the sphaleron 0{x) which 
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is also an extremum of the Euclidean action S'e[0]. We expect this part to yield an 
imaginary contribution to the partition function since we want to describe thermal 
transitions by the sphaleron. Assuming \Zi \ <^ Zq^ we may write 



lm{f^TF) = -lm(lnZ) ^ -Im 



In ( Zo ( 1 



1 



Im(Zi 



(33) 



The evaluation around each of the two classical solutions <^j, z = 0, 1 is done in a 
parametrization similar to but now with spacetime dependent parameter fields 
/(r,a;) = (/(ll)(r,a;),/(^)(r,x)): 



Vl + (/(^)(r,x)r 



/ sin/(ll)(r,x) 
\ cos /'•"^(r, x) 



(34) 



In terms of the parameter fields f{T,x), the classical solutions (pi{x) are given by 
/(II)(t,x) = f,{x), /(^)(r,x) = with /o(x) = 0, /i(x) = f{x) = 4arctan fe^^ 



We add fluctuations v{t,x) = [y'^^^'^{T,x),v^-^\T,x)j to both parameter fields to ob- 
tain the fluctuation ansatz 



{t,x) 



( sin (^fi{x) +t;(ll)(r,a;)) ^ 

v'^-^\t, x) 
\ cos(/i(x) + f (II) (r,a;)) / 



"')|'t 'r\\ \ sin fi{x) ^ 



Vl + (^W(r,x)) 

' sin /j(x) + f *^ll)(r, x) cos T, xjj + yu^"'{T, X 

v^-^\t, x) 

/.(2:) -t;(ll)(r,x) sin/,(x) - | ((t;W(r,x))' + (t;(ll)(T, x))') cos/i(x 



(35) 



cos 



ifi{x) + {if{v))) {t,x) + 0Uv^^\\v^^^ 



/ 

(36) 
(37) 



(One should note that we are considering time dependent fluctuations!). We do not 
attach indices to the fluctuations since it is clear that the fluctuations are different 
for 00 and (fi. 
Here, 



v)) {t,x) 

' t;(ll)(r,x)cos/i(x) - i (^(t;(^)(r,z))% (z;(ll)(r,x))^^ sin/i(a;) ^ 

v^^^ij, x) 

— v*^ll)(r, x) sin/j(x) — i ^t'(-'-)(r, x)^ + (^v'-ll)(r, x)^ j cos/i(x) 



(3^ 
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are the fluctuations around the solution ansatz ipi{x) = [sin fi{x), 0, cos fi{x)), con- 
structed such that they respect the sigma model constraint 0^ = 1 (at least up two 
second order in ifj (v)). They are parametrized in terms of two fluctuation functions 
v{t,x) = (^v^-^\t,x),v^'^^\t,x)^ which represent the two dimensions of a unit sphere. 

We now insert (|36D into the Lagrangian (|I]) and expand it to the mentioned second 
order in the fluctuations which yields the linearized Lagrangian (around the classical 
solution 0i) 



Ci{v,v',v) = Ai 



where 

Ui {v, v) 



^(/;(^))'+^v(cos/.(x)) 



- Ui {v, V 



(39) 



+ 



1 
2 



f-{x) -/iV^3(cos/i(x))sin/i(x 
{f-{x)f + /iV^3(cos fi{x)) cos fi{x 



V^3 (cos /i (s) ) COS /i (x) - V^3^3 (cos /i (x) ) sin^ fi{x) (y 



(40) 



where V^s is the derivative of V with respect to (p^. Defining the linearized conjugate 
momenta 



vr(±) 

^(11) 



dC 
dC 



-- X^k,{x)v^^\ h{x):=l + iy\f:{x)f 



(41) 
(42) 



it is now easy to write down the linearized Hamiltonian in terms of (parameter) fields 
and conjugate momenta n = (vr(||), 7r(_L)^ , i.e. 



7ii{v,tf,'K) = TV ■ V — Ci{v,v ,v) 



(43) 



To insert (^3[) into (p7|), we also have to take care of the change in the integration 
measure. We are performing a substitution 0(r, x) v(t,x). Since we now use 
parameter fields which respect the constraint \<f)\ = 1, we get rid of the delta dis- 
tribution in the integral. From the general theory of distributions, it is well-known 



13 



that for the parametrization cj) = (p{f), 



i=l 



A 



l + (/W(r,a;)^))"^ fd/(ll)(r, x)! fd/(^)(r,x)l (44) 



The change to the fluctuations as integration variables is simply a shift, /^"^(r, x) 
ti*^l')(r, x) = /'■"^(r, x) — fi{x), f^-^\T,x) (— > f^-'-^(r, x) = f^-^\T,x), and we have 



1 + (t;(^)(r,x) 

(_r, xj 



dv^^\T,x)\ [rft;(ll)(r,x) 



(45) 



neglecting corrections of order O (J^v^-^\t,x)"^ ^. We also neglect corrections from 

the substitution n(r, x) i— > vf (r, x) which is a more subtle point. 

Therefore, we are now left with two contributions to the partition function, given 
by the integrals {i = 0, 1) 



dv^^\T,x)] [dw(ll)(r,x) 



d7r(_L)(r, x) 



d7r(j|)(r,x) 



xe 



■ J^'^ dr J°°^ dx (^ — inv+'Hi{v,v' ,tt)^ 



(46) 



The Hamiltonian in ( |46| ) is of standard type ( pOj) with an extra factor to the mass 
M = Ai (which reduces to 1 for the vacuum solution 0o given by fo{x) = 1). Thus 
we may perform the standard momentum integration (150) which yields 



V I ^ki{x)v^^^ (r, x)\ V (r, x) } e 



(47) 



where the Euclidean linearized action is given by (with partial integrations) 

r/3T 



Ai 



dr 

I3t 







dr 



I3t 



dxCe,ii^,v',v) 
1 



dx 



Ux) {v^^)Y + (v^W))' 



+ Ui{v,v) 



Xi I dr I dx v*-"''(t, x) /j"(x) + /i^V03(cos /j(x)) sin /j(x) 
dx«(l')(r,x)^lf [T, Xj 

dxv^^\T,x)gfv^^^ 



/3y roo 

dr 



dr 



(48) 
(49) 

(50) 



14 



For the vacuum we have S'e[<^o] = 0, and 

■1 



2(/'(x))2 + /iV(cos/(a:)) 



dx (51) 



is the Euchdean action of the (static !) sphaleron which obviously equals its (static) 
energy (D times /3r, Se[(pi\ = SsphiPr) = Pt£[0i\ = Pr^sph- 
The part linear in the fluctuations, (^9]) , vanishes if 

/;'(x) -/iV^3(cos/,(x))sin/,(x) = (52) 
/;'(x)-/i^in/,(x) = 0, f,':=^ (53) 

after inserting the potential V{(f)^) = 1 — (f)^. f\ This yields the vacuum (z = 0) and 
the well-known classical static solution {i = 1) (pUj). 

The quadratic contributions (|50D are written in terms of the (full, Euclidean) 
Gaussian fluctuation operators 

= - -Q-^ - - /"^ { V^z (cos fi{x)) COS /i (x) - V^s^A (cos fi (x) ) sin^ /i (x) } 

+ /x2cos/,(x) (54) 



(9r^ dx 

Qf = -{l + ^'(/;(x))n|^-^-{/^V,3(cos/,(x))cos/,(x) + (/;(x))2} 

= -{l + /^^(/;(x))n|^-^ + A.2(3cos/.(x)-2) (55) 

where we inserted V{(j)^) = 1 — 0^. These operators describe fluctuations parallel and 
perpendicular to the solutions 0i. 



Writing the deflnition of operator determinants |28 



P{^(ll)(r,x)}e-^^^/o'"'^-n^^'^''''(-)^ir'^^'''(-) = ^4=^,(56) 



det ( 



v/^.W(r,x)|e-^^r^^n^^^'"'(-)^'"^^^"'(-) = , ^ (57) 

det {g?) 



V 



(the "mass" M = Ai is absorbed in the path integration with measure T>{. . .} de- 



pending on Ai, and we will see why the measure factor Jki{x) in ( [57| ) is necessary 



* Eq. (15^ ) can also be solved for other potentials of the class 1^(0'^) = (1 — '/>^)'°, k € N. For 
k = 2, e.g., one obtains f{x) ~ 2arccot (v^/^a^)- 

Also, the fluctuation equations for these potentials can be obtained easily from the general form 
of the equations given in the text. But already for fc = 2, this yields eigenvalue equations to which 
no explicit solutions are known. 
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for some similar reasons) yields in the one-loop approximation 



det (Qf) det (Qf) 



(5^ 



Inserting these results into (^61) finally leads to the following general formula for the 
transition rate: 



r = - 



Im 



TT 



\ 



det 



(0)' 



(1) 



det 



(59) 



We therefore have to evaluate the determinants of the Euclidean Gaussian fiuctuation 
operators, using appropriate regularisation techniques. 



+ 1 



(60) 



4 Fluctuations around the vacuum 

We first consider the vacuum contribution i- e. z = in all formulae of section 
^. This yields 5*6 [t^o] = and the measure factor ki^{x) = 1, and the fiuctuation 
operators are simply 

with y = fxx, fx > 0. Q^^^ has only a continuous spectrum given by 

(^)' + /i^(l + ^^), keM, (61) 

the corresponding generalised eigenfunctions (which are not square integrable) being 

^n,fc(r, x) = e i/jkdxx) with ijk{y) = e'^^ . 

Decomposing the fiuctuations terms of the eigenfunctions \l/„_fc(r, x) 

(which is commonly called "Fourier-Matsubara decomposition" [^]) thus yields 

1 

/i){.(r,x)}e-^^^/o'^^^r.^^^(-)^<°M.-) = n n lf^V+/^'(i+^') 

keMn=-oo [\ Pt J 

(62) 

V = f v^-^\ This infinite product has to be regularised in a convenient way; usually 
this is done by zeta function techniques [^]. For fixed k, one can perform the 
n = —oo... cxD-product which yields the partition function of an harmonic oscillator 
with frequency u{k) := /xa/I + k"^, 

^(^(fc)):= , ^.^^1(^)^3. = exp{-PTHk';PT,f^)), (63) 

$(A:2;/3r,Ai) := ^VlT^ + (l - e-"^'''^) . (64) 

2 pt 
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The notation /^t,/^) = ^) + ^$2(^^; /^tA*) with 

$i(A;2; ^) = ^yiTP, $2(fc2; a) = In (l - g-"^^) (65) 

will be useful in the following. 

What remains to be done is the product over all k & JR. Taking the logarithm of 

(|6^), we write it as an integral over k weighted with the density of the continuous 

vacuum fluctuation states po (which will be discussed later). The final result for the 

vacuum determinants is thus 

1 1 / f°° \ 

exp [-Pt dk poik)^^; (3t, Ai')j (66) 



det ^det 

5 Parallel fluctuation contributions to the transi- 
tion rate 



decomposition, i.e. 



The determinant of ^|j can also be calculated with a standard Fourier-Matsubara 



+ n\\ (67) 



has eigenfunctions ^'^[^(r, x) = e '^i/jI^\x) with eigenvalues 



27m \ ^ 



(^f)' (68) 



\Pt J 

where (uj^^^ are the eigenvalues of 

Hll^il^ = + /^'{l - 2sech^ (^x)} = {4!^y (69) 

which is just the static fluctuation equation (0). This is a special case of eq. 
( p. 701) discussed in Appendix with [/ = s(s + 1) = 2, s = = 1 and E = 

{uj^^^" — 1 which is known as Poschl- Teller equation It has a normalised zero 



mode iI)q^\x) = y^sech (/ix), (cjq"'') = 0, and a continuum (ui^^^^^^ (k) = /i^(l + 
k'^). Decomposing the parallel fluctuations v*^"-* into eigenstates of G\i^\ v^^^\t,x) = 



1 y AW) (f x) vields 



V 



< v^^^'[T, x)> e 2 Jo J-oo ^ ' ' II ^ ' ' 



m n=— CO 
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where m now runs over the values m = and m = k, k ^ M. Again (assuming zeta 
function regularisation), the n-product yields one harmonic oscillator for each fixed 
m. 

To illustrate this appearence of harmonic oscillators, one can also decompose the 
parallel fluctuations in terms of the eigenfunctions of eq. ( |69D which can be normalized 
with respect to the scalar product 



(6,6)||:=/ Ci{x)U^)dx (71) 

on L'^{1R): (4'n\ipm'^ = Snm (with the usual generalisation for scattering states 
which we assume throughout in the following). This decomposition 

v^^Hr,x)=j:c^lHT)ij^JlHx) (72) 

m 

is related to the Fourier-Matsubara decomposition used above by 

1 CO 

E ee^^ (73) 



T 



n=—oo 



For fixed r, (|72D defines a mapping L'^{1R) £ which is an isometry with respect to 
the L?'{]R) scalar product (,)||. As this scalar product corresponds to the measure 
in the fluctuation path integral (|70D , we may use (^) as substitution in the integral 
without unit Jacobian. This yields 



L — 'J 



W V{cl\r)]e ^ ^ ^ 



= n -m^ (74) 

m 2sinh^^^^ 

Each c^-* -integral thus contributes one harmonic oscillator with frequency lo^ . We 
will use this decomposition technique to evaluate the perpendicular fluctuation in- 
tegral (^) where the Fourier-Matsubara decomposition of the fluctuations is not so 
obvious. Details of the substitution used in (|7^ can therefore be found in section 

In (|70D, the only difference compared with the vacuum calculation is the zero 
eigenvalue m = of TYy. This yields a zero mode factor p9| in the determinant of 



the parallel sphaleron fluctuations which is the partition function of an harmonic 
oscillator with zero frequency, i. e. that of a free mode which has to be treated 
separately. This zero mode factor originates from the translation symmetry of the 
model which is broken by the classical solution (p{x) which is discussed in Appendix 
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^ The result of this (formal) evaluation of the n-product in ( |70| ) for m = is thus 
the replacement 



2 



ntL\\pTj voyj 2sinh5;^^ " ^ V^/^r 

where we introduced a finite space length L to get rid of the volume divergence. 
(Later we will divide the transition rate by L to obtain the transition rate per unit 
volume which is independent of this procedure.) 

The continuous product over G iR is performed in exactly the same way as in 
Section ^ and yields 

n , ■ ^ = exp (-Pt r dkp^^{kMe;(3r,f^)) (76) 

fceiR2smh 

where $(/c^; /Sy/i) is defined in (0) and p\\{k) is the density of the continuum states 
of Till. 

Together with eq. (IBS) we thus have 



det = LA^iexp (^-(3t j^Jk [p\^{k) - p,{k)) ^(A:^ /3t, //)) . (77) 

This "parallel fluctuation contribution" to the transition rate ( |59D is a real quantity; 
we thus expect the perpendicular contribution to be imaginary. 

The integral in ([77|) , containing the difference of the densities of continuum states 
of two Hamiltonians, can be evaluated using standard techniques discussed in refs. 
@, m, H, H, |T], H]. The main formula is 

^ dk (pii(fc) - po(^)) ^k^: f^T,y) = ^j ^ dk $(^2. p)—5\\{k) (78) 

where 5\\ {k) is the phase shift between incoming and outgoing asymptotic plane waves 
in the scattering states of eq. ( |5UD as discussed in Appendix ^, eq. ( |197D : S\\{k) = 
Si{k) = — 2arctan(A;) since s = sy = 1. 

Inserting (^) into the integral of ( [77|) yields for the continuum fluctuation con- 
tributions 

II /"OO fj Ir 1 

dk{p{k)-po{k))<^(e;PT,p) = -^ dk^== + —h\\{PTp) (79) 



271 J-oo y/T+¥ Pt 

The first integral in (^) is logarithmically divergent; it is the sum of zero-point 
energy contributions of the oscillators in the Fourier-Matsubara decomposition (|70D 
and has to be renormalised. We ignore this here. 
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The function /ly is given by 

1 r .u^n....dS^\ik) 



/ill (a) = — / dk^2{k;a 
2tx J-oo 



oo dk 
1 foc In f 1 - e-'^^''i+^ 

The integral is convergent for all a > 0, and 

(/?T/i) ~ - ln(/9T^) > (81) 

in the high temperature limit jSr —>■ {a = /^^/i — > 0) which is shown in Appendix 
0, eq. (iOl). 

Finally, collecting all results, the parallel fluctuation contribution to the transition 
rate (^) reads 



det -iv 1 = LA/jie-^ii^^^'^). (82) 



6 Perpendicular fluctuation contributions to the 
transition rate 

In order to calculate the remaining perpendicular fluctuation contribution to the 
transition rate (|59|) , we have to evaluate the perpendicular fluctuation integral (|57D 
{i = 1) around the sphaleron. Here the measure factor (^T]), 

ki{x) = k{x) := 1 + u'^if'ix)f = 1 + 4:K\ech^ (fix) (83) 

with = fi'^u'^ is the important new ingredient which contains the influence of the 
Skyrme term on the model. The fluctuation integral reads 

with 

g^ = -k{x)—+n^, Tt;^ = -— + /i2|l-6sech2(/ix)} (85) 



where H± is the static Gaussian fluctuation operator (|25|) of the static energy func- 
tional. 

6.1 The measure— modified harmonic oscillator decomposi- 
tion 



To obtain the desired Fourier-Matsubara decomposition in (B4), we solve the follow- 
ing "measure-modified static fluctuation equation" 

, 2 



n^i/j(^\x) = (J^^) k{x)ij^^\x) (86) 
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which is a general Sturm-Liouville problem on L'^{]R) with the metric factor k{x) 
on the right hand side multiplying the eigenvalue |Q. (This equation is different 
from (pSl) . We are no longer considering the sphaleron in the space of static field 
configurations where it was first constructed, cf. Section but embedded it into 
the full euclidean field theory!) From Sturm-Liouville theory, we know that the 
eigenfunctions of (|86|) are orthogonal with respect to the scalar product 



/oo 
-oo 



(87) 



We may thus orthonormalize the eigenfunctions Hxipl^^ = (t^^^) k{.x)i'm'' of (|86D, 

{ip^\ip\'^'') ±_ = Sim and expand the fluctuations v'^-^^ in a Fourier-series with r- 
dependent periodic coefficients: 



(This expansion again contains generalised scattering eigenfunctions in the usual 
sense) . 

The fluctuations are real, so that J2m CmK'^)'^rn\^) = 2m (cm-*) (t) ('0.^-') (a^)- 
Inserting (R3) into the integrand of (E^j yields 



[r, X) 



a) 



m 



X) 



\ m / 

X \-k{xr<i^\r)4'-\x) + [4'-^)" k{x)c\'-\M'-\x] 
;(V'i^))*(x)#Hx)A;(x)(ix], 



m I 



thus 



/3t 



dr 



EE 

m i 



dxv(^)(r,x)^l^^i;(^)(r,x) 



/3t 



^9) 
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^ Jo 



(90) 



This is a sum over harmonic oscillators. In particular, there is no contribution from 
the measure factor k in this sum (it disappears with the help of Sturm-Liouville 
theory) . 

The mapping v^'^^ i— > for fixed r is an isometry L^(IR) i"^ with respect to 
the scalar product containing the measure factor k{x) which also appears in the 
path integral measure, so we may substitute (pHj) into (pi) with unit Jacobian: 



V 



n/^{4^Hr)} 



-(.L^')'|cW(r) 



(91) 



Now it is easy to obtain the Fourier-Matsubara decomposition of the perpendic- 
ular fluctuations. With a Fourier expansion of the coefficients c^-'(t). 



1 oo - . 



(92) 



we have 



y- (±) 



X) 



(93) 



similar to that of the parallel fluctuations. In terms of the c^^, the operator Q^f^ is 
completely diagonalised, 

jyxv^-\r,x)g^V-\r,x) = E^J + (4»^)* ' 04) 

and the substitution v^-^\t,x) ^— > is again an isometry with respect to the path 
integral measure, thus 



V 



{ -^W^^ (r, X) } e-^^ n 



n ^^KJje 



( 2-n-n I 
/3t J 



m n=— oo 



= n n 

with unit Jacobian. 



m n=— oo 



2?™ 



(95) 
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Assuming zeta function regularisation we first perform the n-product to 
obtain 

J I V J ^2 sinh ^^^^ 

This result is also clear from (|9lD , but to controll the necessary renormalisation steps 
(which is not our aim here), one has to return to the Fourier-Matsubara decomposi- 
tion. 



6.2 Solution and discussion of the measure— modified static 
fluctuation equation 



In order to use the harmonic oscillator decompostion (p6D, we have to solve (pGl). 
This Sturm-Liouville equation may be rewritten as (with y = /ix) 



+ 



6 + 4/^2 (u;W)' 



sech^?/ } tlj^^\y) = 0. 



(97) 



which is eq. (|170|) with 



E 
U 



6 + 4i/2 (J^^ 



\ 



[UJ 



(9^ 



1 



It is convenient to discuss the results in terms of := ^ (to^-^A where the constant 



Ao 



does not depend on the Skyrme couphng constant A2, so that it is possible 
to discuss the Skyrme-less limit A2 — * in terms of a;^ which will depend only on 
the one parameter = jj^u"^ = 

There are discrete eigenvalues for all m G Wq with 
T 



'25 

T 



r,2 

m 



m + - ) > 



(100) 



and U (^^{uj^-^^^ ^ = 6 + Ak^uo"^ > which are conditions for the existence of an 

eigenvalue which depend on the eigenvalue itself. If a discrete eigenvalue exists, it is 
given by (|182|) which means it is a solution of 



1 

m H — 
2 



(101) 
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We first prove that there are no solutions of ( |101| ) with m > 2 which satisfy ( |100| ). 
Assume that such a solution (m > 2) exists. If it is positive, tD^ > 0, then 
1 — < 1, but the right hand side of eq. ( p.01| ) may be estimated from below: 



> 



m H — 



[2 - m]2 > 1 (102) 



for m > 2. So has to be negative or zero, cD^ < 0, contrary to the assumption. If 
(D^ < is so negative that 6 + Ak^lj"^ < 0, it cannot be a discrete eigenvalue because 
then t/(co'm) < 0- Otherwise, we may estimate ( |100D from above: 



/25 



+ Ak'^uI - m + - < 



— m + 



2-m < 



(103) 



for m > 2. Therefore, a solution of (|100|) and ( |101| ) cannot exist for m > 2. 

Next, we solve (|101| ) naively by squaring the square root (we have to check whether 
the solutions we obtain really solve (|101|) ) to get a quadratic equation 



A{K'){uy + B{m, n')Cj' + C{m) = 



with 



(104) 



C(m) 



2\2 



(1 + 4k 

— (1 + 4/t^) + 2 ( m+ - ) (I-Ak 



25 



— m 



2\ 2 



25 



+ m + 



(105) 



The general solution of (|104|) is 



-i?(m, k ) ± jA{m, K 



2A(«:2) 



where A(m,/t^) is the discriminant 



A(m, 



4 m + - 
' 2 



25 + 4|29-4(m+- 



+ 64(k 



2\2 



It is easy to check that A(m, k^) is nonnegative if > for m = 0, 1, 2. 
Next, we insert the solutions ( |106| ), i. e. 



~,2 



21 / W {I -4/«2 

4(1 + 4^:2) + l,"" + 2 



1 + 4/?' 



2^2 



± 



(m+ i 

:i+4/€2)2^ 



25 + 4|29-4(m+- 



(106) 



(107) 



k2 + 64(^2)2 (108) 
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into ( |101| ) and ( |100| ) to determine the allowed signs (we already remarked that we 
have to check the solutions). We have 



8 (m + i) ± J 25 + 4 (^29 - 4 (m + if) ^ 54(^2)2 



4 '"'^ 2(l+4/€2) 



(109) 



The modulus is always positive for the upper (plus) sign, and one can easily check 
that this is a solution of eq. ( |101|) which satisfies ( |100|) for all > 0. But for the 



lower (minus) sign, we have to distinguish two cases: For < < y|- j^-^ — , cu,^ 



does not solve ( |101| ). Otherwise, for > ^ 1 — ^ jt solves this equation, but 

the condition (|100|) is not fulfilled. We therefore have to reject oj'^ _ as solution to 
( |101|) under the condition ( |100|) . So finally, there are three discrete eigenvalues of 
(13) given by 



) - u;,{K } - —-^ - (i + 4«:2)2 

(110) 

^iV) = ^^ = (111) 



, , (4^^)' 4K2-f + f J25 + 16«:2 + 64(/t2)2 

4(/.2)=4('^') = ^^ = ' y.^, ^>0. (112) 

/i^ (1 + 4k^)^ 



There is a zero eigenvalue ('^1''''') = 0. In Appendix 0, we show that this zero mode 
is due to the U{1) rotation symmetry of the Euclidean action ( pi]) which is broken by 
the classical sphaleron solution 0i{x) (analogous to the translation zero mode which 
can be found in the spectrum of the parallel fluctuations). This effect is independent 
of k"^, so the zero mode exists for all .2. 

It is interesting to analyse the .^-dependence of the other two discrete eigenvalues 
which is shown in Fig. |^. For 0, we have cJq(.^ = 0) = —3 which is the Skyrme- 
less result (A2 ^0). In the same limit, we have a)|(.^ = 0) = 1 which explains why 
there are only two discrete eigenvalues in the original model [|^: Without a Skyrme 
term, the third eigenvalue is the lower bound of the continuum. In the limit — > cx), 
all three eigenvalues tend to zero. This can be interpreted as the Ai — limit of the 
model in which there is no sphaleron solution (and therefore no stability spectrum). 

The corresponding eigenfunctions are given by ( |184| ) were s now depends on 

. 2 



(u;^)) (cf. |9|). Using (^), we have 



8m.2 _ 1 + /25 + 4 f 29 - 4 (^m + J .2 + 64(/t2)2 
-U-^^)- ^ ^7^^^;;^^ . (113) 



25 
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Figure 3: Discrete eigenvalues uj'^{k'^) of 



For m = 1, this yields si{k^) 
the eigenfunctions is thus 



2 which does not depend on . The explicit form of 



(y) 



N^'^^sechy tanh y 



1 



2S2(«:2 



(114) 
(115) 

sech^y I (116) 



where Nj;'^^ are appropriate normalisation constants. Like the zero eigenvalue, the 



zero mode ( [L15| ) is independent of 
In the limit 



K 



and agrees with that given in the literature 



2g 



0, we have So(0) = 2, so ipo^\y; k'^ = 0) = sech^y which again 
is Mottola-Wipf 's result. The —>■ limit of 1^2'^^ is even more interesting, because 
there is no such third discrete eigenvalue in the Skyrme-less theory. In fact, we have 
52(0) = 2 and therefore 



lim lim il^2^\y; 



1 



0) L\IR). 



;ii7) 



But since S2(k^) > 2 for > and 'i(^2^\y; n 



exponentially, '^2'''^ is a 



square integrable eigenfunction (and really a discrete eigenvalue) for > 0. 

The continuous spectrum starts from E > 0, i. e. u'^{k) > /i^. We denote the 
continuous eigenvalues by uj'^{k) = /i^(l + A;^), k E IR. The asymptotic behaviour of 
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the scattering states belonging to is discussed in Appendix ^ (|193| - |195| ), Now 

the parameter s is not fixed but depends on k (and k^): 



s = s±{k; K^) = J— + 4/t2(l + A;2) - - e iR, (118) 

and the coefficient of the incoming wave (|195|) is 

2\ r(— zfc)r(l — ik) 

^ ' " ^ " T{-ik - s{k- K-^))V{-ik + s{k- K^) + 1) ■ ^^^^^ 

6.3 Evaluation of the perpendicular fluctuation determinant 



With the above resuhs we perform the discrete part of the product (|96D . From m = 2, 
we simply get a harmonic oscillator factor 

(120) 



2sinh(^^ 



)/3tM 



Next, we consider m = 0. For sufficiently large T, there is only one negative factor 
in the product (^) which means that has exactly one negative eigenvalue. Then 
we may interpret the negative mode as giving rise to an imaginary factor in Zi, 
according to the prescription 

1 11 , , 

121 



2sinh( "-("P^^^ ) 2 2^sh.^ 2 



where the additional factor | arises from the distortion of the non-Gaussian contour 



over half of its range p5| . 

Finally, m = in ( pSj ) is the zero mode contribution of a)f = which is discussed 
in Appendix 0, yielding ( |169| ) 



L— , := 4a/a^,/-^a/i + -H? (122) 

2sinh^i^ ^ iy3/?^^V ^5 > 

(since this zero mode results from the f/(2)-rotation symmetry of the Lagrangian 
broken by the sphaleron, it does not lead to a space volume divergence like the 
translation zero mode factor (|75D; the parameter describing this symmetry takes 
values in [0,27r)); cf. ( |T55|) ). 

The product over the continuum u^'> = uj^-^\k), k G IR in ( PBj ) is rewritten as 
an integral over $(A;2; (3t,1^) defined in (|6^) weighted by the density p_l(A;; k^) of the 
scattering states of eq. (P7|), 

n , = i-^T r dkp^ik;K')^e;PT,t^)) (123) 

fc6iR2smh V -^-^ / 
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We are now able to write down the contribution of the perpendicular fluctuations 
to the transition rate. From (||) with z = 1 and (|96| Jl20|JMl , [l22| , |T23D we have 



det 



(1) 



2 



8zsin(^^-(fi^) sinh(^ 

xexp(^-/3T J^^dk {p^{k;K^)-po{k))^k^;(3T,fi^)) (124) 

The integral is performed using the analogue of formula ([78|), 

dk U{k;K^)-poik))<l>{k';(3T,^i) = 7r dk^e;PT,f^):jrS±{k;f^') 

-oo ^ ' Zn J-oo CLK 

(125) 

where d±{k] k^) is again the phase shift between incoming and outgoing plane waves 
in the scattering states of 



^ 9, r(—ik)r(l — ik) 

5^{k-K^) = arg— — ^ / / J 126 

r(— ^A; — s_l(A;; K^))T[—ik + s<^[k; k^) + 1) 

Using the decomposition ^{k'^] /3t, p) = $i(A:^; /i) + ■g^$2(^^, /^t/^) (pSf), it is easy to 
show that the integral 

/oo rj 
dk^,{k^-p)—5^{k-K^) (127) 
-oo dk 

is again a divergent zero point energy contribution which has to be renormalized. 
With $2(fc^, /^tm); we define the function (a = (3tp) 

h^a, := ^ r In (l - ^-"^) ^^^4^^^ (128) 
27r j-oo ^ ' dk 

which is discussed in Appendix p. Its high temperature behaviour for /t^ > is given 
by (gg), 

/^±(/3t/^;«:')~^^4^7^>0 for/?T/i-.0. (129) 

This linear growth with temperature results from the large k behaviour of the phase 
shift ( pTlD , 

- 5^{k]K^)^ m{K^)k + b{K^) for A; ^ oo. (130) 

In the Skyrme-less limit {k^ = ^ sx_{k] = 0) = 2), this phase shift tends towards 
a constant, and the high temperature leading term of the continuum fluctuation 
contributions is again logarithmic ( p06| ), 

/ix(a; = 0) 2 ln(/3T/i) > for (3tP — ^ 0. (131) 
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Finally, the perpendicular fluctuation contribution to the transition rate (|59| ) is 



det 



\ \gi^ I 82 sin (h^MMdf) sinh 



(132) 



7 The thermal transition rate due to the sphale- 
ron: Skyrme— less limit and range of validity of 
the calculation 

So far we have evaluated explicitly the one-loop fluctuation contributions to the 
perturbative transition rate (|59|) . The classical (zeroth order) part of this relation, 



Im 



TT 



\ 



det 



(0)' 



(1) 



det 



(133) 



is given by the exponential of the classical action of the sphaleron (^Tj), i.e. 



r 



i(/(x))2 + ^V(cos/(a:)) 



dx = S/^T/iAi. (134) 



Inserting now (|8^ , 132 ) and ( |134|) into ( |133| ) (where the damping factor |^| for large 
Ai — which corresponds to the weak coupling limit of the original 0(3)-cr model — 
may be approximated by the modulus of the negative eigenvalue, \g\ = |(jj_(k^)|), the 
transition rate per unit volume given by the sphaleron is 



r _ |u;_(/t^)|/i 
L ~ n 



exp |— 8/9t/^Ai — h{/3TfJ.', «;^)| 

(135) 



with the continuum fluctuation contributions described by the function 



(136) 



The prefactors of the exponential originate from the transition rate damping con- 
stant (flrst factor), the zero modes (numerator of the second factor) and the nonzero 
discrete modes of the fluctuation spectrum (denominator of the second factor). 

The Skyrme-less limit of the model is given by — 0. In this limit, there is no 
discrete eigenmode with positive eigenvalue a;+(K^), the contribution ^ 



being absorbed in the continuum fluctuation contribution. 



exp{-/t(/3r^;/t^)} ^2^0 r , 2 n^l 

2dnh( -^^-;^^-^ ^ -""^P|-^(^^^'^ =0)|. 



(137) 
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With = 0)) = —3, the Skyrme-less transition rate is given by 



^""'^ ^ exp{-8/5T/iAi - /i(/?t^;k2 = 0)} (138) 



L 



_ 2Ai// 1 

k2=o 



7r/3T sin 



which is exactly the result of Mottola and Wipf (eq. (5.10) in their paper), using our 
notation (Ai ^ p-, o;^ yu^). 

Next, we discuss the temperature range in which the sphaleron calculation is valid. 
We will see that this range is influenced by the Skyrme term. Since we performed a 
one-loop quantum approximation around the classical sphaleron solution using the 
idea that the imaginary part of the partition function describes transition rates, this 
range of validity is restricted by two technical requirements: 

1. There has to be exactly one negative eigenvalue of the full Gaussian fluctuation 
operator. This leads to a lower bound. 

2. The one loop fluctuation contributions must not be greater than the zeroth 
approximation given by the classical Boltzmann term. This leads to an upper 
bound. 



The negative eigenvalue of the static fluctuation operator is (|110|) , i.e 



and the corresponding eigenvalues of the full Gaussian fluctuation operator are given 
by 

'27rn\ / , 2n\2 2 



+ [uj4k')} (140) 

n = always yields a negative eigenvalue, the eigenvalues for n > thus have to be 
nonnegative. Inserting n = 1, one obtains from condition |I| 

kBT> '^-(/)'^ =: ksT^'\^,n'). (141) 

Without Skyrme term, this lower bound is T^^\fi,K^ = 0) = -^/x. With Skyrme 
term, it decreases with increasing Skyrme coupling (since |u;_(/t2)| decreases with 
increasing k). But for these lower temperatures, the exponential suppression due to 
the Boltzmann factor increases. 

To analyse the second condition, one has to compare the classical action of the 
sphaleron, 5*6 [<^i] = S/^y/iAi with the quantum effects of the linear fluctuations around 



it, given by the continuum fluctuation function h in the exponent of (|135|) and the 
zero and discrete eigenmode factors in that formula. In fact, it is sufficient to compare 
5'e[<^i] with hi^pTli] n'^), because for high temperatures (small Pt), the discrete mode 



30 



prefactors of ( |135| ) yield only logarithmic contributions in the exponent, and h 
itself has at least a logarithmic /^t/^ — ^ behaviour. The upper temperature bound 
is thus given by the condition 

8/?T//Ai > hiPrl^; «;^) (142) 

in which we ignore all numerical prefactors to obtain a rough estimate of the upper 
temperature bound T^^^/i, Ai; k^). 

In the Skyrme-less case = 0, the high temperature behaviour of h is given by 
(Pip, ( |131D , i.e. /i(/3T/i, = 0) ~ — 31n(/3r/i), 0. For small arguments, any 



negative power function grows faster than the negative logarithm, lima^o l^^"'*'* = 
for all 7 > 0. This means that for any 7 > 0, there must exist a number a* (7) such 
that for all a < a* (7), we have — ln(a) < a~'^. Inserting this into ( |142|) without 
numerical factors yields 

keT < //Ap (143) 

The upper temperature limit therefore increases with Ai: High temperatures are al- 
lowed for large Ai which is the small coupling limit of the original 0(3)-cr model. 
This means we consider small /^tA* = a in the approximation of the logarithm 
— ln(a) < a~'^, therefore small 7 7^ are allowed in this approximation. Thus, 
the upper temperature bound in the "small coupling limit" Ai ^ 00 is given by 

kBT < Ai/i =: kBT^^\Xifi; = 0). (144) 

Next, including the Skyrme term changes the high temperature behaviour of h 
drastically. In view of (|129|) we have 

/i(/3T/i; > 0) ~ '"'^^'''^ } for (3tI^^0. (145) 
6 PrfJ^ 

Inserting this into ( p.42| ) and again ignoring numerical prefactors immediateley yields 



ksT « =: A;bTW(Ai, ^; > 0) (146) 

y m(K^j 

For fixed > 0, this upper bound is considerably smaller than the Skyrme-less 
upper bound (|143|) in the "small coupling limit". One main effect of the Skyrme 
term in the sphaleron calculation of the thermal transition rate is thus to reduce the 
range of validity of the calculation. In the Skyrme-less limit, T'^+^(Ai,/i; > 0) in 
eq. ( |145| ) increases as m(/t^) decreases for — which leads back to the Skyrme-less 



upper temperature bound condition (|143|) . 



8 Discussion and conclusions 

After the detailed discussion of the sphaleron and the fluctuations around it in the 
Skyrmed reduced 0{3)-a model we enquire about the physical implications of the 
results. 
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First of all, from a very naive point of view, one might conclude that eq. ( |135| ) 
describes thermal transitions over the barrier separating the degenerate vacua which 
are exponentially suppressed at high temperatures. This would destroy the most im- 
portant property of the sphaleron induced transition — namely that of not leading 
to exponential suppression — compared to that induced by the (vacuum) instanton 
which is always exponentially suppressed. But eq. ( p.35| ) does not demand this in- 
terpretation, because in the temperature region where this exponential suppression 
becomes dominant, this formula is no longer valid for the reasons discussed in the 
last section. 

At very high temperatures, T > T^~^\ the sphaleron calculation thus breaks down. 
In realistic models with spontaneous symmetry breaking (and particularly in the 
electroweak model), this is cured by the restoration of symmetry above Tj-est ~ 
(which is indicated in the diagrams of Fig. |l| and Fig. |^). At these temperatures, 
the fluctuations become very large and the system becomes subject to disorder. It is 
clear that in this temperature range the transition rate becomes unsuppressed and 
there is no need to worry about the sphaleron which ceases to exist because it needs 
the broken symmetry. In 0(3)-cr models (with or without Skyrme term), where 
there is no Higgs-mechanism of spontaneous symmetry breaking, one has to add 
a symmetry breaking potential iJ,V{(j)^) by hand. One can show that the explicitly 
broken symmetry of the 0{3)-a model is not restored at high temperatures, so that 
for very high temperatures T > T^~^\ both the Skyrme-less and the Skyrmed reduced 
0(3)-cr model cease to be physically meaningful theories. The only difference in the 
very high temperature behaviour is therefore that T^"*"^ decreases with increasing 
Skyrme coupling k^, but this is what one expects on physical grounds, because at 
high temperatures, the Skyrme term (having the highest number of derivatives in the 
Lagrangian) dominates the theory with large fluctuation effects which would lead to 
symmetry restoration also at lower temperatures in a realistic theory. 

Finally, we consider the situation in the lower temperature region around T^~\ 
As already remarked in the introduction, the Skyrme-less limit of the model (i.e. the 
original reduced 0(3)-cr model) is a type II theory with dominant transition processes 
as indicated in Fig. This changes completely for > since the Skyrme term 
enables the existence of localised instantons with Euclidean energy (^ |T5[] 

2So = 8nXiw{K), w{k) = 1 + 0{k) (147) 

(the function w{k) has to be evaluated numerically since the vacuum instanton is not 
known explicitly ||2l[]) and corresponding periodic instantons with < O.In the 

limit — *■ 0, there is thus a discontinuous change in the type of instantons involved 
(whereas the action of the vacuum instanton is continuous for ^ 0, the action of 
the constrained instanton in the Skyrme-less theory can be estimated to be 2>S'o = 
SttAi 13). In the case of the sphaleron of the Skyrmed model, there is a smooth 



crossover to the Skyrme-less situation. The sphaleron is only quantitatively changed 
by the Skyrme term, whereas the instanton changes qualitatively. In particular, the 
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Euclidean action S^pl of the sphaleron at the lower temperature limit \ 



depends smoothly on the Skyrme coupling k"^. This means that the type II condition 

2So<si;l ^ ro>T(-) (149) 

which is true for = remains true for > due to the smooth K^-dependence 
of both actions, whereas the type of instantons changes discontinuously from those 
in type II theories to those in type I theories. The Skyrmed reduced 0{3)-a model 
is therefore neither of type I nor of type II at least for a weak Skyrme coupling |jl5| . 
Whether it may become a type I theory with increasing Skyrme coupling depends 
on the function w{k,) which would have to satisfy the condition w{k)\uj-{k'^)\ > 2 for 
some K^. But this is very unlikely to happen, because as shown in Fig. H, |u;_(k^)| 
decreases rapidly to as increases. 
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Appendix 

A Zero modes of the 0(3) model 

The spectra of the operators Tiy , 7i_L each contain one zero eigenvalue which appears 
in the evaluation of the fluctuation integrals (|^J37|) and yields divergences in (|7D|,[ 



In this Appendix, we discuss the origin of these "zero modes" p9 |. 

For this reason we analyse the symmetries of the model (|l]) which are broken 
by the classical solution 0i{x) = (sin /(x), 0, cos /(x)). We consider the following 
symmetries of the model ([^): 

1. space translations: 

T[(/), a](r, x) = 0(r, a; + a) 

= (p{T,x) + ^{t,x) ■ a + 0{a^) for a^O (150) 
ox 

(There is also a time translation symmetry, but this is not broken by a static 
solution.) 



33 



2. internal S'0(2)-rotations: 



R[0, a](r, x) 




T,x) ■ a + 0{a ) for a — 

(151) 



The classical solution 0i{x) = (sin /(x), 0, cos /(x)) breaks this invariance because 
T[(^i, a]{T, x) = (pi{x + a) ^ (fi{x), R[(^i, a](r, x) 7^ (^i(x). (152) 



The moduli space 

is therefore topologically A4 = IR^ x and one can show that ^[(pi,Q](x) and 



of the parameters a, a describing finite energy solutions 

9T 

da 



^[(^1, 0](x) are zero modes of the static Gaussian fluctuation operator 
by expanding 



£ 



and £ 



9^i + ^[c^i,0] 



(153) 



around (pi. The problem here is that one has to be very careful with the constraint 
0^ = 1, and it is better to consider the broken symmetries in the parametrization 



(^(/). Therefore we first express the symmetries ( p.5CI| ), ( |151| ) in terms of the param- 
eter fields. This is done by expanding the ansatz (|^) in small quantities 5f around 
f = (/^"''i Z*^"'"'') and then comparing the result with the infinitesimal form of the 
symmetries given by ( [1501) , (|151|) , also written in terms of the parameter fields. The 
result is 



1. for space translations: 



T/[/, a](r, x) = /(r, x) + ^(r, z) ■ a + (9(a^) for a 



(154) 



2. for internal S'(9(2)-rotations: 



-/W(r,x) cos/(ll)(r,x) 
(l + (/W(r,x))2) sin/(ll)(r,x) 



for a^O (155) 
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In terms of the parametrization /, the classical solution ip is given by f{x) = fo{x) = 
{f{x), 0). This yields the following two zero modes of the Gaussian fluctuation oper- 
ator in /-parametrization, 



[/o,0] = ; , fix) = 2/isech(/ix) (156) 



da 

[fo,0] = ( ^ ) , sin/(a;) = — 2sech(yua;) tanh(yua;) (157) 




da 

Obviously, they are proportional to the zero eigenvalue eigenfunctions of the fluctu- 
ation equations (p9|j86|) , normalised with respect to the appropriate scalar products 

= i/fsech(/xx), = 1 (158) 

'^i^\x) = ]J ^^^_^^^^^ sech{fix) tanh(/xx), = 1 (159) 

Next, we need a convenient way to take these zero modes into account in the 



Fourier-Matsubara evaluation of the fluctuation integrals (|70| , |95|) . Usually, on uses 
the method of collective coordinates which means that one treats the parameters 
which describe the breakdown of a global symmetry of the theory in a static solution 
(a and a in the case of the 0(3) model) as additional dynamical degrees of freedom, 
i.e. one considers 

${t,x) = R[T[(^,a(T)],a(r)](r,x) (160) 

and thus adds two redundant degrees of freedom a, a. One can show that the re- 
sulting theory is a gauge theory, i. e. the Euclidean action of the theory in terms 
of the new "fields" Se[v,a,a] {v denoting again the parameter fields of the fluctua- 
tions) is invariant under some (time-)local transformations of a, a. Therefore, the 
usual techniques for path-integral quantization of gauge theories (Faddeev-Popov 
37| or BRST |3^) can be applied to evaluate the partition function Zi around cpi. 



A convenient gauge-fixing condition is to exclude fluctuations parallel to the zero 
modes. This procedure is well-known and discussed in detail in the hterature, even 
in the particular context of 0(3) models, using classical Faddeev-Popov or more 



modern BRST and BV methods Therefore, we do not discuss the details. 



Here, we use a simpler semiclassical method [20] to compute the zero mode con- 



tributions to the partition function, i.e. to treat them in fluctuation integrals ([70i95). 
In the harmonic oscillator decomposition of the Gaussian fluctuation operators G^^\ 

G±\ each zero mode contributes one harmonic oscillator mode with zero frequency, 
i. e. a free mode. The classical contribution of a free mode to the partition function 
is 

^e~^^'4 = ^ f dg. (161) 
27r V2^J ^ ^ ' 
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The "coordinate" q of this free mode has to be identified with the zero mode coef- 
ficients Cg"'', ci'^ in the harmonic oscillator decompositions of the parallel and per- 
pendicular fluctuations, (|7^ and (pH]), respectively. Therefore, we consider the fluc- 
tuations in direction of the zero modes. Let prg be the projector from the space of 
fluctuations v{t, x) on the subspace of zero mode fluctuations. Since the zero modes 
are those which result from fluctuations in the collective coordinates, expanding the 



(162) 



parameter field version of ( |160| ) 

/(r, x) = R/[Ty[/;, a{T)la{r)]{T., x) 
around the classical static solution /o = R/[Tj[/o, 0], 0] (small a, a) yields 



R/[T/[/o,0],0](r,x) 



da 



[7^,0], 



(r, x) ■ a(T) + 



rift 

^[T,[/;,0],0](r,x) 



with 



= fo{x) + 
= fo{x) + \/Xipro{v{T,x)} 

prQ{v{T,x)} = -7^[io,0](r,z) ■ a(r) + —-^[fo,0]{T,x) ■ ^(t) 



"^^^ [fo, 0](r, x) ■ a(r) + ^[fo, 0](t, x) ■ a(r) 



da 



da 



(163) 



(164) 



da da 
(The factor Ai appears because an overall factor Ai is extracted from the fluctuation 
expansion (^8|-^)). 

But from the decompositions ( |72D and ([88|), the parallel and perpendicular com 
ponents of the projection operator pr^ are given by prg"'' = (•j'i/'o 



ipQ^^ and 



(-L) 

pro 



Vi^^) so that with ([T5^,|l5 



Aipro{v(T,x)} = JXi 



cS\t)^?{x) 



x] 



a{r)f'{x) 
Q;(r) sin f{x) 



Thus we have 
dq 



I 



dq 



Ai 
A^ 



dCn 



dci 



AiV(/',/')||Aa 



with 



(/',/')ii = V 



Ai Ac\ 



(sin /, sin/)_ 



Ai y (sin/, sin/)_LAa 
10 + 16k^ 4 



(165) 

(166) 
(167) 

(168) 



15 /i 

Aa, Aa are the values the parameters describing the broken symmetry can take: 
For rotations, it is clear that Aa = 2n. To avoid a volume divergence due to the 
translation parameter a G iR, we restrict the model to a finite space volume such 
that Aa = L. Therefore the zero mode factors are 



= 4VAi 



71" 



3/3T/i 



'1 



(169) 
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B The Schrodinger equation with sech —potential 



All attempts to evaluate the fluctuation factors (^,^) lead to the eigenvalue differ- 
ential equation 

tp" + {E + Usech'^y)'ilj = (170) 

with ip = ip{y) which we discuss in this section |^2||. There are three different cases 
to be distinguished: 

1. f/ > 0: then it is possible to have 

(a) < 0, i. e. bound states (discrete spectrum) 

(b) > 0, i. e. scattering states (continuous spectrum) 

2. f/ < 0, then E > Q and only scattering states (continuous spectrum) exist. 
We first set ^ = tanhy to get (now ip = ip{C,))'- 

E 



d_ 



%p = Q 



(171) 



Assuming ?7 > 0, we first compute the discrete spectrum for which i? < is valid. 
We may thus define 

e := \/^, e > (172) 

and 

1 /l 
^ = -2 



'- + f/>0 



U = s{s + l), S-- 
to obtain from ( |171| ) the differential equation of generalized Legendre functions 

^2 



(173) 



d 



+ 



s(s + l) - 



1-e 



ip = 



(174) 



Next we set ^p{^) = (1 — ^'^)^^w{^) which yields 



d^ 



2 2(e + l)^ + [s(s + l)-e-e> = 



(175) 



Finally, changing the argument to x = |(1 ~ 0) ( P-75| ) transforms to 



d'^w 



dw 



Xix -l)-r^ + [{e + l)- (2e + 2)x]— - {e - s){e + s + l)w = 



dx 



dx 



(176) 



which is a hypergeometric differential equation. 

The general form of the hypergeometric differential equation is 

ziz - 1)^ + [(c - (a + 6 + l)z]^ -ab<j) = 
dz^ dz 



(177) 
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Solutions are given by the hypergeometric functions 



F{a,b,c; z) = J2 

k=0 



{a)k{b)k 



with 



T{a + k) 
Via) 



a{a + 1) ■ . . . ■ (a + — 1), (a)c 



1. 



(178) 



(179) 



The series ( [178|) converges for |2;| < 1, 2; G (T. It also converges for 2; = 1 if c—a—b > 0. 

Comparing ( |176D with the general form ( |177D , we identify a = e — s,b = e + s + l 
and c = e + 1, so the solution to eq. ( |176|) is given by 



wix) = F{e- s,e + s + l,e + l;x)■ 



:m 



Inserting all the substitutions, the solution to the original equation ( p.70| ) may be 
written as 



ijj^y) = (1 - tanh^y)' F (e - s,e + s + l,e + 1; ^(1 - tanhy)) . 



;i81) 



We want to compute bound states (discrete eigenvalues for E < 0). Hence the 
solutions (|180|) have to be square integrable which means that at least (as all functions 

are continuous) i/j{y) — ?^ 0. In terms of the x-variable, which is related to the 
original variable y hj x = |(1 ~ tanhy), this means that w{x) has to be bounded 

as x ^' — ^ and ^ — >° 1. We always have F{a, 6, c; 0) = 1, so the limit y ^ cxd is no 
problem. But since c — a — b = —e < 0, F{a, b, c; 1) is not bounded for general values 
of a, b, c. We thus have to force F(e — s, e + s + 1, e + 1; 1) < 00 which is achieved by 
making the infinite sum ( |178D finite: (m)^ = for — m G JNq and k > m. This yields 
the quantization condition e — s = —m G JNq: m labels the eigenvalues of (|1 74|) , 
^m = s — m. There are no eigenvalues for any m G INq in view of the condition ( |172D 
e > which means m < s. The discrete spectrum of eq. ( |17CI| ) (case is therefore 
given by 




with 



1 
2 



me INn and Trjim) > 



Tu(m) 



4 2 



and the corresponding eigenfunctions are (with 1 — tanh^ y = sech^y) 
ipm{y) = sech.^^~"^\y)F (—m, 2s — m + 1, s — m + 1; -(1 — tanhy) 



;i82) 



;i83) 



;i84) 
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Next, we assume U < 0. Then, there is only a positive, continuous spectrum 
fcase 0) with E > 0, and we substitute 



k := ±VE =^ {±iky = -E (185) 
U = s{s + l), s = ~ + ^^ + U e(U. (186) 

Performing the same steps as above, the scattering states of eq. ( |170| ) for ?7 < are 
given by 



ik 



'^±,k{y) ~ ij- ~ tanh yj ^ F \^±ik — s, ±ik + s + 1, ±ik + 1; -(1 — tanhy) 

(187) 

and the continuous spectrum is 

{e = k'^\k e R} . (188) 

It is obvious that these calculations are also true for U > 0, E > (then, s is real), 
so that (|187|) also describes the scattering states of the U > potential which has 



bound and scattering eigenstates (case fT^). 

Next, we analyse the asymptotic behaviour of the scattering states (|187|) . It will 
be sufficient to consider the scattering state i/j^^kiu)- The prefactor in ( p.87| ) behaves 
like 

\ - tanh^ y)~~ ~ 4^e^''= for x — > ±oo (189) 



. ik 



(we ignore the common factor 4~ in the following.) To deal with the hypergeometric 
function, we use ( [178| ), i. e. F{a,b,c; z) = 1 + 0{z) for z ^ Q. This works well for 
y +00, because then |(1 — tanh?y) — > 0. For y — > — oo, we use the following 
identity: 

F{a,b,c;z) = r(c)r(c - a - b) ^ ^ ^ ^ ^ ^ _ ^. ^ _ 
1 (c — a)L [c — b) 

+ - ^y-'-'Fic -a,c-b,c+l-a-b;l-z) 

r(a)r(6) 

(190) 

which yields, applied to ( |187|) 



F (^—ik — s, —ik + s + 1, —ik + 1; -(1 — tanhy)^ 



T{-ik + l)T{ik) ^ /. , ., , , 1,, 

^ ^ ^ -F i-ik - s,-ik + s + l,-ik + + tsinhy) 



+ 



r(-s)r(i + s) V , ^ . ^ . , . 
T(-ik + i)T(-ik) fi^^ , 

— — ^ ' ; '- - 1 + tanh2/ 

V{-ik - s)V{-ik + S + 1) 

xF (^s + 1, -s, ik + 1; ^(1 + tanhy)^ . (191) 
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Inserting 

^(1 ± tanhy) ~ e^'^y for y — ^ ^oo (192) 

and putting everything together, we get the following asymptotic behaviour of the 
scattering states: 

^-^k{y) A{k)e-'''y + B{k)e'''y (193) 

^_,,(2/) '^rt^ e^^y. (194) 

with 

r(— s)r(l + s) T[—ik — s)T[—ik + s+l) 

Since k"^ > 0, the poles of B{k) are excluded. A{k) vanishes for s G JN. This is the 
well-known case of the reflection free Rosen-Morse potential. 

Eqs. (|193| , |194|) show that i/j-^k describes incoming particles from the left; there 



are no reflected particles to the right (one can check that the second independent 
solution describes incoming particles from the right). The incoming flux is 

normalised to unity. Thus, all important information about the scattering solutions 
resides in the phase shift S defined by the transmission coefficient B(k) = e^^: 

r r T(-ik)T(l-ik) , , 

T{—tk — sjT{—ik + s + 1) 

For s G IV, there is a simple formula for 6s{k) ||3l| : 

6s{k) = -2j2 arctan . (197) 

C (High) temperature behaviour of the functions 

h\\ and h i 



In this Appendix, we analyse the small a = f^TjJ' behaviour of the functions /ly (^ 
h± ( |128| ) given by the integrals 



- r In (l - e-^) ^^^dk (198) 

with (IT96D 



,,,, T(-ik)T(l-ik) , , 

i y—%k — s)V y—ik + s + 1) 

bg is known for arbitrary values of s, but we will only use it for integer s E JN (cf. 
formula ([T97|)). 



40 



For the parallel fluctuation continuum, we have s\\ = 1 and ^^l;^ 



hu (a) 



1 j"^ 

TV J~oo 



In f 1 - e-'''^^ 



dk. 



so that 



(200) 



1 + P 

This integral is convergent for a > 0. For high temperatures we have a = 0. 
Using the approximation 



In (l - e-'^^^ 
the substitution Q = ak yields 

h\\{a) ^ - 

which means 



In (aVl + = - In (a^ + {ak)' 



(201) 
(202) 

/i|| (/^T/i) ~ -ln(/3r^) > for [3t^ — > 0. (203) 

In the case of perpendicular fluctuations, the behaviour of the phase shift is much 
more complicated. We have 5s^{k) = 6±{k; k^) because (cf. eq. ( p,18| )) 



a r"^ In (a^ + O 
2^ loo a? + C 



ln(2a) 



s_L = s_l(A;; 



'2^ 1 
-+4/.2(i + P)_-. 



(204) 



The Skyrme-less limit {k^ = 0) reduces this to s_l(A;;k^ = 0) = 2 and 



dk 



~i+k^ ~ 4H^' '^^^ same approximations as in the parallel case which 

yield 



h±{a; 



K 



0) = / In (1 - e-" 

rin(a^ + C 
- (ln(2a) + ln(3a)) 



C,=ak 



1 



1 + A;2 4 + A;2 
2 

+ 



dk 



and thus 



/i^(a; = 0) ~ -2 ln(/3T/i) > for prf^ 



0. 



(205) 



(206) 



Including the Skyrme term {k^ > 0) changes this behaviour. First, in order to 
restrict the following discussion to positive k, we perform an integration by parts in 
( p.98| ) and use the symmetry of the integrand to obtain 



hi (a: K, 



1 

TT Jo 



-S±{k; 



d_ r, 
dk 



In (l - e-" 



dk. 



(207) 



For flxed k, s±{k\K^) and thus —5_i{k]K^) increase with go hj_{a; k!^) increases 
with increasing Skyrme coupling for flxed a. Now the phase shift (|196|) , i.e. 

Sj_(k; K^) = arg v / v / 



F [-zk + \ - y^f + 4/€2(l + F)) F [-ik + \ + y^f + 4/€2(l + P) 

(20 
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can be evaluated numerically. The main result is that for fixed k^, one can always 
find positive real numbers a = a^K!^) > 1, c = c(k^) > such that 

• for A; G (0,a), \S±{k; k^)] is bounded by c < cx3 

• for k > a, —6±{k; k^) increases linearly. 

It is this linear asymptotic behaviour of 5±{k\ k^) = 0{k) for k ^ oo which deter- 
mines the leading a-term in the asymptotic a —>■ 0-expansion of (|207D . Here, we 
see one major infiuence of the Skyrme term on the original 0{3)-a model, because 
without Skyrme-term, we have S±{k; = 0) = O (^-p^ for k oo. 

Besides a numerical analysis, one can also visualize the linear asymptotic be- 
haviour of 6±{k] K^) by some simple analytic arguments. Consider the monotonically 
increasing series (fcj) of real positive numbers for which s±{ki; /t^) G IN. For large fcj, 
we may estimate s±{ki; k^) ~ [nki] where the Gauss bracket [] denotes the integer 
part of a real number. With 



arctan 



K 



fki' 

< arctan — 
\ n . 
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n G {1, . . . , [i^ki]}, 



one can use ( |1971 ) to obtain 

1 

arctan 



1 



[Kki] < -6±{ki; K^)<- - {[Kh] + 1) 



(209) 



(210) 



Thus, the series (— 5±(A;i; K^))jejv increases linearly for i ^ oo. It should be pos- 
sible to use the general formula for 6±{ki; h^) to show that this is also true for the 
interpolating values of k. 

We may therefore use the approximation 



< 



-5^{k; 



c(/€2) 

m{K^)k + b{K^) 



for < k < a 
for k > a 



(211) 

is difficult to 
> 27r for 



with m{K^) = 0{k) for — > oo. The small K^-behaviour of m[K, 
estimate, but it is clear from the above discussion that m{K^) — > 0, b{K^^ 

^ 0; the Skyrme-less limit thus changes the asymptotic k behaviour of the phase 
shift S±{k; K^) from a linear decrease to a constant value. 

With ( |211|) , we estimate the integral (P07|) : 



< 



1 

TT Jo 



-6±{k;^'] 
1 



d 

dk 



TC 



In (l - e-'^^^ 



— e 



dk 



dk. (212) 



The first integral is easy: 



TT 



dk 



In 1 - e 



dk 



C[K 



TT 



In 



1 - e 



Infv^ 



a' 



1 - e-"" 
+ 0{a) 



for a 0.(213) 
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The second integral splits up into two parts, 



OO fj 

k— 

TV Ja dk 



In (l 



dk + 



°° d 



TT Ja dk 



ln(l 



dk (214) 



with 



p ^ [in (l - e-"^'^)] dk = - In [l - e-''v^i+"')) = C(ln(a)) for 



la dk 

In order to evaluate the remaining contribution, we insert 



a ^ 0. 
(215) 



0< Aln(l-e-"^ 



k 



a 



< a 



dk V ^ ; Vl + A;2 e«^^i+^ - 1 e"'^' - 1 

into the first integral of ( pi4| ). Using the dilogarithmic function 

Ht) 



for A; > 1 (216) 



dilog(x) : = 



1 l-t 



-dt 



we obtain 



''dk 



ln(l 



dk ~ 



ak 



1 CdC 



-dk 



a Jaa — 1 

^ + -dilog (e"-^) + ^a'a 
ba a 



1 



71 



a + 0(a) for a 

6a 



0. 



(217) 



(218) 



Obviously, this is the leading term in the asymptotic a —>■ expansion of hj_{a] k^) 
Thus 



7rm(K 



2^ I 



for /^T/i — 0. 



(219) 



6 /^T/i 

In the Skyrme-less limit, this term vanishes {m{K^ = 0) = 0), and the leading 
behaviour of h±{PTfJ'] 0) is given by ( P06|) . 
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